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1.0 INTRODUCTION 


This report contains a description and user's guide of the 
conputer program developed at VPI&SU under the NASA contract NAS9-15091. 
The major goal of tills cratract was the development of a FORTRAN com- 
puter program to simulate the dynamics of the NASA-DELCO Electro- 
mechanical Actuator for aerospace applications. The program developed 
under this contract had to meet the following criteria: 

1. Adequate accounting of the effects of the stator 
phase currents on the permanent magnets of the 
rotor. 

2. Reasonably accurate simulation of the voltage 
and current waveforms present In the power con- 
ditioner network during the motoring, regenera- 
tive braking, and plugging nodes of operation. 

Such accuracy requires; 

a) Accounting for the major component non- 
llnearlcles such as the coupling or 
chopper Inductor saturation and the diode 
and transistor nonlinearities, 

and b) Inclusion of the machine nonlinearities 
(see Item 1 above) . 

3. Reasonably accurate simulation of the overall EMA 
servo loop during a step response. This requires 
the Inclusion of all major nonllnearltles such as 
machine velocity saturation, machine torque limiting, 
etc. 


These requirements have all been met by the computer model 
developed under this contract. Tlie derivation of this model Is developed 
in detail In the first volume of this report. 
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2.0 MODEL DEVELOPMENT 

ITiie four channel EM actiiator shown in Figure (2-i) is usually 
operated with two machines or channels active while the ether machines 
are braked. Therefore, one needs to consider only two channels at a 
tine in the modeling process. These two channels do not operate inde- 
pendently of each other because of the velocity correction controller 
which keeps both machines nmnlng at approximately the same speed. 

Since both channels operate under very similar conditions it 
was decided to model the one and two machine modes based upon the 
single channel mechanization diagram shown in Figure (2-2) . Under 
this assumption the velocity correction signal shown in this diagram 
would always be zero. Figure (2-3) shows essentially the same 
mechanization diagram with the addition of the blocks representing 
the rotating masses. The differentiation between. one and two channel 
operation is provided by the ficticious gear ratio NX. When only 
one machine is operating NX equals 2 since the machine has to rotate 
through twice the angular displacement when compared with the normal 
two machine mode. With two machines operating NX is reduced to 1. 

Closer inspection of this diagram reveals that there are 
essentially three basic components to this one and two channel model. 
They are: 

1. The power conditioner and machine 

2. The velocity correction loop 

3. The position correction-mechanical loop 

Notice that these components have been listed in order of ascending 
time constants. A schematic of the power conditioner and machine 


network (item 1) ia given in Figure C2-4) . The electrical model is 
shown in Figure (2“’5). Figure (2-6) displays the various components 
of the EMA vhJ.ch make up the system of rotating masses. 

The models for each of the EMA components are derived in the 
standard state space form symbolized by 
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gO!<Ea CONDITIONER AND MACHINE NETWORK STATE EQUATIONS 

Ihtt state equacloos of the power conditioner and oachlne network can 
be written as follows: 


A X + B u 
— n -n — n —n 


where 


^n 


F • G 
-n -n 


2n 


• H 

-n -n 


These matrices and vectors are defined as follows: 
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The coefficients Cll, C12, ••• C34 in matrix [G^] relate the branch 
currents IBIS, IB19, and IB20 to the state variables XSV as follows 
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3.0 PROGRAM DESCRIPTION 

A program vas wrlcten In FORTRAN to solve ntsserlcally the 
state equations of the three subsystems presented In the previous 
chapter. Since the time constants for these three componencs vary 
widely. It was decided to integrate each one of these loops Indepen- 
dently. For example, the mechanical-position loop which varies the 
slowest TiTould be Integrated first. The velocity loop which has the 
next smallest time step would then be Integrated NIirr2 times for 
each, of the NINTl time steps of the mechanical-position loop. The 
velocity loop integrations use linear Interpolation on the old and 
new state vectors of the mechanical loop in order to generate the 
Inputs. PE and R7EL. The same is true for the power conditioner and 
machine network. In this case, there are NINT3 integrations per 
velocity loop integration. This scheme eliminates having t.o 
Integrate all 14 state variables at a rate determined by the smallest 
time step. The relationship between the three separate integration 
loops is given in the following example: 

1) Let NINTl - 3 

NINT2 « 2 
NINT3 - 2 
and TSNET - 10 ps 

2) Then TSVLP - NINT3*TSNET - 20 ps 


^ _ 


and TSMLP - NI1IT2*TSVLP - 40 ps 
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iple, the relationship between the 
loops is given below: 
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In order to allow some flexibility In the computer simulations 
It was decided to have three different inodes of program execution. 

The first mode, PQIODE ■ 1, considers the power conditioner and machine 
network only. The other two modes include the entire servo loop 
but with a simple machine model, for PGMODE " 2 and a detailed model 
for PGSiODE » 3. A flow chart of the overall logic flow of this pro- 
gram la given in Figure C3-1) • 




23 


PROGRAM MOEE 1 (PGMODE - 1) 

ThJLs mode Is characterized by the following: 

1. Fourth order machine and power conditioner model. 

2. Constant machine speed Cl*e., there are no mechanical loop 
integrations) 

3. Constant torque (current) command (l.e., there are no velocity 
loop integrations) 

Typical Applications 

This mode is useful when it is desired to examine the behavior 
of the power conditioner and machine network at a specified speed, 
torque command and mode of operation (i*e., motoring, regenerative 
braking, and plugging) . Total simulation time is generally below 
10 ms (EMA time) which corresponds to 1000 integrations with a time 
step, TSMET, of 10 ys. 


PROGRAM MODE 2 (PGMODE = 2) 


This mode is characterized by the following: 

1. Simple machine and power conditioner model consisting of: 

a) machine torque, TM = KT * IM, and 

b) a machine speed characteristic which is clamped at 
zVMSAT (typically 9500 rpm) 

2. Fifth order velocity loop model 

3. Fifth order mechanical loop model 

Typical Applications 

This mode is used to obtain the step response of the EMA with 
.one or two machines operating. The use of the simple machine model 
greatly reduces the total execution time and at the same time gives 
excellent results when compared to test. 
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PROGRAM MODE 3 (PGMODE » 3) 

This mode is characterized by the following: 

1. Fourth order machine and power conditioner model. 

2. Fifth order velocity loop model. 

3. Fifth order mechanical loop model. 


Typical Applications 

This mode is designed to obtain step responses of the EMA with 
either one or two machines operational. The fourth order power con- 
ditioner and machine model generates more accurate machine torques 
which serve as inputs to the mechanical-position loop model. Typical 
simulation times CEMA, reference) are generally of the order of 
40,000 integrations of the network state equations. This large number 
of integrations is necessary because of the rapid response of the 
power conditioner and machine model when compared with the velocity 
and mechanical loop integrations. 
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4.0 RUN PREPARATION 


This chapter describes the steps in the preparation of the 
input data cards for the EHA simulation program. The input structure 
has been designed to allow the execution of several cases per run. 
This input data structure can be represented plctorlally as follows: 


NCASES 


CARD #1 




ICASE = 1 
RUN DATA ' 


■' S| 

ICASE = 2 
iRUN DATA 

O 

O 

o 

ICASE = NCASES 
RUN DATA 


SET n 


SET #2 


SET # NCASES 


INPUT 

DATA 

DECK 


The total number of cases per run is given on the first data card 
and is represented by the integer variable NCASES in the program 
Itself. The card sets which follow contain the input data for the 
cases specified by NCASES. In the main program each individual 


case is identified by the variable ICASE. The required variables 
and their input formats are the topics of the remainder of this 
chapter. 
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DESCRIPTION OF INPUT DATA 


1. NCASES - Number of simulations or cases per run 


PGMODE - Program Mode 


P(^ODE 


1. Machine Model Only - 4th Order 

2. Complete EMA Model - lOth Order 

(Simple Machine Model) 

3. Complete EMA Model - 14th Order 

(Detailed Machine Model) 


3- lOMDDE - Printing control 


lOMODE 


0 Do not print results of integrations 

1 Print results of every integration 

* O 

0 
0 

N Print results of every Nth integration 


4. IDISK - Disk storage control 


IDISK 


0 Do not load integration data on disk 

1 Load specified data for each integration on disk 

O 

O 

o 

N Load specified data for every Nth integration on disk 


5. 

6 . 
7. 


J 0 Do not plot 
1 1 Plot data stored on disk 


IPLT - Line printer plot control 
IPLT 

NPLTS - Number of X-Y Plots (Excluding the plot of the logic signals) 
IPX - Location in INDEX of the X variable for the X-Y Plots 
IDX - NPLTS + 1 


8. NPAGES - Number of pages per line printer plot 

9. NVSTOR - Number of X-Y values stored on^'disk 

NVSTOR = NPLTS + 1 

10. INDEX - Vector which contains the ID numbers of the X-Y variables 
to be stored on disk. The ID numbers of the Y variables 
come first. The ID number of the X variable is placed 
immediately following the last Y variable. (See list 
of ID codes later in this report) 








11. NIRT3 - Number of Network Integrations: 


PQIODE » 1 - Total Number of Network Integrations 
PGMODE - 2 - NINT3 » 0 

PG^K)DE “ 3 - Total Number of Network Integrations 
per Velocity Loop Integration. 


12. NTERMS - Control of the transition matrix calculation for the 
Network model 


NTERMS = 


0 Use the error tolerance ETNET to control 
PHI and THETA 

>1 Number of terms for PHI and THETA 


13. TSNET - Network integration time step (seconds) 

5 ps £ TSNET £ 20 ps 

14. ETNET - Maximum allowed difference in the last two terms of 

PHI for the PC-Machine Network 

15. SHIFT - Commutation Switch Control (mech. rad.) 

16. ITOL - Machine Line Current Tolerance (amps) 

17. ^ - Machine Torque Constant (nt - m/amp) 

18. SAMPLE - Chopper Inductance Control 


SAMPLE 


0 Linear Case: Do not update inductance 

of the chopper inductance 
0 < SAMPLE <1 e.g. SAMPLE = 0.5: Update 

Chopper inductance in 5% 
increments 


19. XSV(l) • • • XSV(4) - Power Conditioner and Machine Network State 

Variables 

20. ^ - Supply Battery Voltage (volts) 

21. RDl* • • RDM - On resistances of diodes D1 - DM (ohms) 

22. ROFF - Equivalent resistance of a cut off diode or transistor (ohms) 

23. RQl- ♦ • ROM - On resistances of transistors Q1 - QM (ohms) 

24. LI, L3, L4 - Leakage inductances of phases a, b, and c respectively. 

(Henries) 


25. RIO - Battery Resistance (Ohms) 


26. 


27. 


! 8 . 


29. 

30. 

31. 


32. 


33. 


34 . 


35. 


36. 

37. 


38. 


39. 


40. 


41. 


42. 


43. 


44. 


R3 - Chopper Inductor Resistance (ohms) 

L2 - Chopper Inductor Inductance (henries) 

(Maximum Value) 

Cl - Capacitance of Filter Capacitor 
(farads) 

NBCHS - Number of Network Branches 
NTOIGS - Number of Network Twigs 
NLINKS - Number of Links 
NSTWIG - Number of Twig Voltage Sources 
NCTWI G - Number of Twig Capacitors 
NRTWIG - Number of Twig Resistors 
NLTWIG - Number of Twig Inductors 
NCLINK - Number of Link Capacitors 
NRLINK - Number of Link Resistors 
NLLINK - Number of Link Inductors 

NDPTS - Number of data points in the curve of the chopper 
Incremental inductance 

XI(1) • • XI (NDPTS) - X-values of the chopper incremental inductance 

curve (amps) 

FXI (1) • • • FXI (NDPTS) - Y-values of the chopper incremental inductance 

curve (henries) 

NINT2 - Number of velocity loop integrations per mechanical- 
position loop integration 

IF PGMODE = 1, NINT2 = 0 

NTERM2 - Control of the transition matrix calculation for the 
velocity loop model 


NTERM2 


0 Use the error tolerance ETVLP 
to control PHIV and THETAV 
^1 Number of terms for PHIV and THETAV 


TSVLP - Time step of the velocity loop (seconds) 


TSVLP - NINT3 * TSNET 


45. 


ETVLP - Maximum allowed difference in the last two terms of 
PHIV for the velocity loop model 

46. TAUl - Position error amplifier time constant (seconds) 

47. TAU2 - Velocity error FB amplifier time constant (seconds) 

48. ^ - Velocity error amplifier gain 

49. KRL - Current error integrator gain 

50. - Velocity FB amplifier gain 

(amps /rotor radian/sec) 

51. VHI - Maximum (+) rotor velocity during plugging (radians/sec) 

52. VLO - Maximum (-) rotor velocity during plugging (radians/sec) 

53. FlYl - Maximum plugging current divided by 5 

54. F1Y2 - liaximum non-plugging current divided by 5 

55. F2Y - Current rate limiter Y-coordinate (.machiuc:^ amps/5) 

56. F2X - Current rate limiter X-coordinate (machine amps/5) 

57. :\V(1) - ICMD/5 

58. XV(2) - ICMDl/5 

IF PGMODE = 1 Set XV(2) equal to the desired machine 
line current divided by 5 

59. XV(3) - Velocity FB-Loop Output 

60. XV(4) - ICMDL/5 

XV(5) - Current Rate Limiter Output 
(machine amps * .0005277/sec) 

62. NMACH - Number of machines active (1 or 2) 

63. NINTl - Number of integrations of the mechanical-position loop 

64. NTERMl - Control of the transition matrix calculation for the 

mechanical-position loop model 

f 0 


NTERMl 


>1 


Use the error tolerance ETMLP 
to control PHIM and THETAM 
Number of terms for PHIM and THETAM 
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65. T5MLP - Time step of che mechanical-position loop, (seconds) 

66. ETMLP - Maximum allowed difference in the last two terms of PHIM 

for the mechanical-position loops 

67. TAU3 - Position error amplifier time constant (seconds) 

68. ^ - Deflection command (flap degrees) 

69. VMSAT - Rotor velocity/ limit (rotor radians/sec) 

Applies to PGMODE = 2 only 

2 

70. - Machine and reflected gear inertia (kg - m“) 

2 

71. JF - Flap inertia (kg - m“) 

72. ^ - Flap viscous damping coefficient (nt - m/flap rad/sec) 

73. XACT - Actuator mount stiffness (nt - m/EMA output rad) 

74. ^ - Flap Stiffness Coefficient (nt - m/flap rad) 

75. - Position FB-Loop Gain (deg/ rad) 

76. KERR - Position error amplifier gain 

(machine amps/ 5 /degree flap) 

77. ^ - Gear ratio 

78. ^ - Gear ratio 

79. ^ - Gear ratio 

80. XM(1) - Rotor velocity (rad/sec) 

81. XM(2) - EMA Reaction Torque (nt - m) 

82. :gl(3) - Flap Velocity (rad/sec) 

83. XM(4) - Flap Position (radians) 

84. XM(5) - Amplified Position Error Signal (machine amps/5) 

85. XLAB(l)* • • XLAB(IO) - Contains the x-axis label for plots 

(40 characters max) 

86. YLA3(1)« » ‘ YLAB(IO) - Contains the y-axis label for plots 

(40 characters max) 

87. TITLE(1)« • » TITLE(IO) - Contains the x-y plot title 

(40 characters max) 
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INPUT DATA FORMAT 


The input data is read from 80 column computer cards. The formats 
for each variable are designated as follows: 


VAR *♦ 


no 


■Variable Name 


° specifies the format 
code (110 in this case) 

Starting column of the 10 column field 
corresponding to format 110 


The placement of integer, real single precision, real double precision 
and character data is described by the examples given below. 

Integer data Input format: 110 

e.g. IVAR » 10 and 


no 


Col # 


1 2 3 4 5 , 6 , 7 , 8 ,9. , 10 


DATA CARD: 


n-n 


m. 


*RIGHT HAND JUSTIFY ALL INTEGER DATA 
Real Single Precision Data Input Format E14.7 
e.g. RVi\R = 10. 


and 


E14.7 


Form 1: 


Form 2: 


12 3 

4 

5 

6 7 8 

Q 

> 


11 12 

13 

14 • • • 

80 

ll 1 0 1 • 

rn 

Z] 

TTZ 

r: 

i 




1 1 


d 

The 10. 

can 

be 

placed anywhere 

within 

the 

14 spaces 


12 3 


5. „ 

6 7 3 


10 , 

11 12 

L2 

14 • ■ • 

80 


n 



Cl 

□ 


E 

1 1 • • • 

i 1 


Real Double Precision Data 


e.g. DVAR = 10. DO 


Input Format D14.7 

(i.e. 10) and |l j D14.7 


1 

2 

3 


5 _ 


7 , 





19 13 

14 

• * • ftO 

Ci] 

□ 

□ 

□ 


□ 

d 

□ 

□ 

□ 

□ 

1° 

1 ^ 

-inJ 


*Notice the exponent D1 must be right hand justified. 














FORMATS FOR THE INPUT DATA CARDS 

1. NCASES 

iif iioj 

2. PGMODE lOMODE IDISK IPLT 
[j] IIP I 111 no 1211 IIP l31l ifot 

3. (SKIP IF DISK » 0) 

NP1 .TS IPX NPAGES NVSTOR 
nr I107ll{ IIP i2li IIP bli IIP "1 


(SKIP IF DISK -= 0) 


INDEX(l) ••• 

m no ncn no" 


1211 no l3ir no 1411 no isii no ieil ii 


5. NINT3 NTERMS 

111 no iiii no' 


6. TSNET 

HI wrrr 


7. SAMPLE 
111 ei7 


ETNET 

fill eiT! 


SHIFT ITOL KT 

j29| E14.7 M E14.7 j57| zVTT 



8. XSV(l) 


XSV(2) 

iiM^ 


XSV(3) XSV(4) 

1291 D14.7 I43t D: 


9. E4 

111 E1477 


RD2 RD3 RD4 RD5 

1151 S14.7 1291 E14.7 l43i E14.7 l57i E14.7I 


11. RD6 RPR RDB 

HI E14.7 IlM S14.7 1291 E14.7" 


12. RQI RQ2 RQ3 R04 

HI ElO nil E14.7 1291 E14.7 S43l E14.7 


El^.7 
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STORAGE OF PLOT DATA ON DISK 

All data which is to be plotted is first stored on a temporary 
disk file inorder to reduce the core requirements of the program. 

This data is loaded onto aisk using the following formats: 


1. 5E14.7 

(Real data) 

2. 80L1 

(Logical Variables) 


The unit number of this temporary disk file at VFI is 8. If this 
unit number is different at the USER'S installation then the following 
changes must be made: 

In' subroutines 

PDATA 

EMAPLT 

BOLPLT 

all statements of the form 

REWINDS and 
WRITE (8, ) 

must be changed to the unit number valid at that 
installation. 

For example if the unit number is 9 then these become: 

REWIND9 

WRITE (9, FORMAT #) 

This data when read from disk by the plotting routines is loaded into 
plot vectors which presently have a dimension of 1001 words. This 
corresponds to the data produced by 1000 integrations (i.e. Initial 
Conditions + 1000 integrations). If more than a 1000 integrations 
are contemplated then some of the data must be skipped. This is 
accomplished by the input parameter IDISK defined earlier in this 


chapter. 
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IDENTIFICATION NUI'SERS FOR VARIABLES WHICH MAY BE STORED ON THE 
TEMPORARY DISK FILE 

(THESE ID NUMBERS ARE LOADED INTO THE VECTOR INDEX) 

In addition to these ID numbers, each variable also has a restriction 
code number defined as follows: 

MEANING 

No restrictions 
Not available for PGMODE =» 1 
Not available for PGMODE = 2 
Not available for PGMODE » 3 
Not presently available 

Do not plot or store on disk (meaningless) 

The ID numbers and restriction codes for each variable are listed 
sequentially starting on the following page. 

The loading of the variables identified in INDEX on to the temporary 
disk file is accomplished by SUBROUTINE- PDATA. 


RESTRICTION CODE 

0 

1 

2 

3 

4 

5 






45 


VARIABLE NAME 

ID NUMBER 

RESTRICTION CODE 

LOCATION 

CIA 

1 

2 

COMMON/BLKl/ 

CIB 

2 

2 

COMMON /BLBa/ 

CIC 

3 

2 

COMMON/BLKl/ 

VA 

4 

2 

COMMDN/BLKl/ 

VB 

5 

2 

COMMON/BLKl/ 

VC 

6 

2 

COMMON/BLKl/ 

VAB 

7 

2 

COMMON/BLKl/ 

\'BC 

8 

2 

COMMON/BLKl/ 

VGA 

9 

2 

COMI^ON/BLKl/ 

PA 

io 

2 

COMMON/BLKl/ 

PB 

11 

2 

COMMON/BLKl/ 

PC 

12 

2 

COMMON/BLKl/ 

PCORE 

13 

0 

COMMON/BLKl/ 

PUL 

14 

0 

COMMON/BLKl/ 

PTSM 

15 

2 

COMMON/BLKl/ 

PSO 

16 

2 

COMMON/BLKl/ 

PEM 

17 

2 

COMMON/BLKl/ 

PGM 

18 

0 

COMMON/BLKl/ 

PNM 

19 

0 

COMMON /BLiCl/ 

TEM 

20 

2 

COMMON/BLKl/ 

TM ■ 

21 

2 

COMMON/BLKl/ 

RRANG 

22 

2 

COMMON/BLK6/ 

RANG 

23 

0 

COMMON/BLK6/ 

RVEL 

24 

0 

COMMON/BLK6/ 

RACEL 

25 

4 

COMMON /BLK6/ 

ICMDL 

26 

0 

COMMON/BLK7/ 

ICMDl 

27 

0 

COMMON/BLK7/ 

IMC 

28 

0 

COMMDN/BLK7/ 

IM 

29 

0 

COMMON/BLK7/ 
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VARIABLE NAME 

ID NUMBER 

RESTRICTION CODE 

LOCATION 

BOC(l) 

30 

2 

COMMON/BLK12/ 

BOC(2) 

31 

2 

COMMON/BLK12/ 

BOC(3) 

32 

2 

C0MM0N/BLK12/ 

BGAP(l) 

33 

2 

C0M1'!0N/BLK12/ 

BGAPC2) 

34 

2 

COMMON/BLIQ2/ 

BGAP(3) 

35 

2 

COMMON /BLK12/ 

FS(1) 

36 

2 

COMMON /BLKl 2/ 

FS(2) 

37 

2 

COMMON/BLK12/ 

FS(3) 

38 

2 

COMMON/BLK12/ 

EM(1) 

39 

2 

COMMON/BLK12/ 

FM(2) 

40 

2 

COMMON/BLK12/ 

FM(3) 

41 

2 

COMMON/BLK12/ 

FG(1) 

42 

2 

COMMON /BLK12/ 

FG(2) 

43 

2 

COMMON/BLK12/ 

FG(3) 

44 

2 

COMMON/BLK12/ 

HM(1) 

45 

2 

COMMON/BLK12/ 

HM(2) 

46 

2 

COMMON/BLBC12/ 

HM(3) 

47 

2 

C0MM0N/3LK12/ 

HG(1) 

48 

2 

C0MM0N/BLK12/ 

HG(2) 

49 

2 

COMMON/BLK12/ 

HG(3) 

50 

2 

COMMON/BLK12/ 

LMV(l) 

51 

2 

COMI^ON/BLK12/ 

LMV(2) 

52 

2 

COMMON/BLK12/ 

LMV(3) 

53 

2 

COMMON/BLK12/ 

LGV(l) 

54 

2 

COMMON/BLK12/ 

LGV(2) 

55 

2 

C0MM0N/BLK12/ 

LGV(3) 

56 

2 

C0MM0N/BLK12/ 
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s 


VARIABLE NAME 

ID NUMBER 

RESTRICTION CODE 

LOCATION 

QLOSS (1) 

57 

2 

COMMON /3LK13/ 

(2) 

58 

2 

COMMON/3LK13/ 

(3) 

59 

2 

COMMON/BLK13/ 

(4) 

60 

2 

COMMON/BLK13/ 

(5) 

61 

2 

COMMON/BLK13/ 

(6) 

62 

2 

COMMON/BLK13/ 

(7) 

63 

2 

COMMON/BLK13/ 

(3) 

64 

2 

COMMON /BLKl 3/ 

DLOSS (1) 

. 55 

2 

COMMON/BLK13/ 

(2) 

66 

2 

COMMON /3LK13/ 

(3) 

67 

2 

C01-1M0N/BLK13/ 

(4) 

68 

2 

COMI'ION/BLKIS/ 

(5) 

69 

2 

COMMON/BLiaS/ 

(6) 

70 

2 

COMMON/BLK13/ 

(7) 

71 

2 

COMMON /BLK13/ 

(8) 

72 

2 

COMMON/BLK13/ 

(9) 

73 

2 

COMMON/BLia3/ 

PELOSS 

74 

2 

COMMON/BLK13/ 

FANG 

75 

1 

COMMON/BLK28/ 

PE 

76 

1 

COMMON/BLK28/ 

VE 

77 

1 

COMMON/BLK28/ 

TACT 

78 

1 

COMMON/BLK28/ 


TIME 


79 


0 


COMMON/BLK27/ 
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VARIABLE TIME 

ID NUMBER 

RESTRICTION CODE 

LOCATION 

BCHCUR (1) 

201 

2 

COMMON/BLK2/ 

(2) 

202 

2 

COM^X)N/BLK2/ 

(3) 

203 

2 

COMMON /BLK2/ 

(4) 

204 

2 

COMMON/BLK2/ 

(5) 

205 

2 

COMMON/BLK2/ 

(6) 

206 

2 

COMMON /BLK2/ 

(7) 

207 

2 

C0MM0N/BLK2/ 

(8) 

208 

2 

COMMON/ BLK2/ 

(9) 

209 

2 

C0M140N/BLK2/ 

(10) 

210 

2 

C0MM0N/BLK2/ 

(11) 

211 

2 

C0MM0N/3LK2/ 

(12) 

212 

2 

C0MM0N/3LK2/ 

(13) 

213 

2 

C0MM0N/BLK2/ 

(14) 

214 

2 

C0MM0N/BLK2/ 

(15) 

215 

2 

C0MM0N/BLK2/ 

(16). 

216 

2 

C0MM0N/BLK2/ 

(17) 

217 

2 

C0MM0N/BLK2/ 

(18) 

218 

2 

C0MM0N/?.LK2/ 

(19) 

219 

2 

COMMON /BLK2/ 

(20) 

220 

2 

C0HM0N/BLK2/ 

(21) 

221 

2 

C0MM0N/BLK2/ 

(22) 

222 

2 

C0MM0N/BLK2/ 

(23) 

223 

2 

C0MM0N/BLK2/ 

BCHEMF (1) 

224 

2 

C0MM0N/BLK2/ 

(2) 

225 

2 

COMMON /BLK2/ 

(3) 

226 

2 

C0MM0N/BLK2/ 

(4) 

227 

2 

C0MM0N/BLK2/ 

(5) 

223 

2 

C0MM0N/BLK2/ 

(6) 

229 

2 

C0MM0N/BLK2/ 

(7) 

230 

2 

C0MM0N/BLK2/ 

(8) 

231 

2 

C0MM0N/BLK2/ 


VARIABLE NAME 

ID NUMBER 

RESTRICTION CODE 

LOCATION 

BCHEMF (9) 

232 

2 

C0MM0N/BLK2/ 

(10) 

233 

2 

C0MI>I0N/BLK2/ 

(11) 

234 

2 

C0M10N/BLK2/ 

(12) 

235 

2 

common/blb:2/ 

(13) 

236 

2 

COMMON/ BLK2/ 

(14) 

237 

2 

C0M1'I0N/BLK2/ 

(15) 

238 

2 

COMMON/ BLK2/ 

(16) 

239 

2 

C0MM0N/BLK2/ 

(17) 

240 

2 

COMMON/ BLK2/ 

(18) 

241 

2 

COMMON/ BLK2/ 

(19) 

242 

2 

COI4MON/BLK2/ 

(20) 

243 

2 

COMMON/ BLK2/ 

(21) . 

244 

2 

COMMON/ BLK2/ 

(22) 

245 

2 

C0M^^0N/BLK2/ 

(23) 

246 

2 

COMMON/ 3LK2/ 

BCHPOW (1) 

247 

2 

COMMON/ BLK2/ 

(2) 

. 248 

2 

COMMON/ BLK2/ 

(3) 

249 

2 

COMMON/ BLK2/ 

(4) 

250 

2 

COMMON/ BLK2/ 

(5) 

251 . 

2 

C0MM0N/BLK2/ 

(6) 

252 

2 

COMMON/ BLK2/ 

(7) 

253 

2 

C0MM0N/BLK2/ 

(8) 

254 

2 

C0MM0N/BLK2/ 

(9) 

255 

2 

COMMON/ BLK2/ 

(10) 

256 

2 

C0MM0N/BLK2/ 

(11) 

257 

2 

C0MM0N/BLK2/ 

(12) 

258 

2 

C0MM0N/BLK2/ 

(13) 

259 

2 

C0MM0N/BLK2/ 

(14) 

260 

2 

COMMON/ BLK2/ 

(15) 

261 

2 

C0MM0N/BLK2/ 

(16) 

262 

2 

C0MM0N/BLK2/ 


VARIABLE NAME 


ID NUMBER 


RESTRICTION CODE 


LOCATION 


BCHPOW(17) 

263 

2 

C0t^M0N/BLK2/ 

(18) 

264 

2 

COMMON/ BLK2/ 

(19) 

265 

2 

COMMON/ BLK2/ 

(20) 

266 

2 

COMMON/ BLK2/ 

(21) 

267 

2 

C0M^[0N/BLK2/ 

(22) 

263 

2 

COMMON/ BLK2/ 

(23) 

269 

2 

C0MM0N/BLK2/ 

XSV(l) 

270 

2 

C0MM0N/BLK15/ 

XSV(2) 

271 

2 

COMMON/ BLK15/ 

:tsv(3) 

272 

2 

COMMON/ BLK15/ 

XSV(4) 

273 

2 

C0M1^0N/3LK15/ 

U(l) 

274 

2 

C0MM0N/BLK15/ 

U(2) 

275 

2 

COMMON/ BLK15/ 

U(3) 

276 

2 

COMMON/ BLIC15/ 

U(4) 

277 

2 

COMMON/ BLK15/ 

XSVDOT(l) 

278 

2 

COMMON/ BLK15/ 

XS\T)0T(2) 

279 

2 

COMMON/ BLK15/ 

XSVD0T(3) 

280 

2 

C0MM0N/BLK15/ 

XSVD0T(4) 

281 

2 

COMMON/ BLK15/ 


XV (1) 

282 

■ * 

1 


C0MM0N/BLK17/ 

(2) 

283 

1 


COMMON/ ELKl 7/ 

(3) 

284 

1 


COMMON/ BLKl 7/ 

(4) 

285 

1 

DO NOT 
USE THESE 

C0MM0N/BLK17/ 

(5) 

286 

1 

• UNLESS 

C0MM0N/BLia7/ 

UV(1) 

287 

1 

NINT2*1 

AND 

C0MM0N/BLIQ7/ 

(2) 

288 

1 

NINT3-1 

C0MM0N/BLK17/ 

(3) 

289 

1 


C0MM0N/BLK17/ 

(4) 

290 

1 


C0MM0N/BLK17/ 

(5) 

291 

1 

J 


C0MM0N/BLK17/ 




VARIABLE NAME 

ID NUMBER RESTRICTION CODE 

LOCATION 

XVOLD(l) 

292 

5 

COMMON/ 3LK1 7/ 

(2) 

293 

5 

COMMON/ BLKl 7/ 

(3) 

294 

5 

COMMON/ BLKl 7/ 

(4) 

295 

5 

COMMON/ BLKl 7/ 

(5) 

296 

5 

COMMCN/3LX17/ 

UVOLD(l) 

297 

4 

COMMON/ BLKl 7/ 

(2) 

, 298 

4 

C0MM0N/BLK17/ 

(3) 

299 

4 

COMMON/ BLK17/ 

(4) 

300 

4 

COMMON/ BLK17/ 

(5) 

301 

4 

COMMON/ BLKl 7/ 

XVI (1) 

302 

1 

COMMON/ BLKl 7/ 

(2) 

303 

1 

C01'1M0N/BLK17/ 

(3) 

304 

1 

C0MMON/BLK17/ 

(4) 

305 

1 

COMMON/ BLKl 7/ 

(5) 

306 

1 

COMMON/BLK17/ 

UVI(l) 

307 

4 

COMMON/ BLKl 7/ 

(2) 

308 

4 

COMMON/ BLKl 7/ 

(3) 

309 

4 

COMMON/ BLKl 7/ 

(^) 

310 

4 

COMMON/BLK17/ 

(5) 

311 

4 

COMMON/BLK17/ 

XM(1) 

312 ■ 

1 

C0MMON/BLK20/ 

(2) 

313 

1 

COMMON/BLK20/ 

(3) 

314 

1 

COMMON/3LK20/ 

(4) 

315 

DO NOT , 

USE THESE 

COMMON/BLK20/ 

(5) 

316 

UNLESS 1 

» 

COMMON/BLK20/ 

0M(1) 

317 

NINT2-1 , 

AND 

COMMON/BLK20/ 

(2) 

318 

NINT3-1 1 

COMMON/BLK20/ 

(3) 

319 

1 

COMMON/BLK20/ 

(4) 

320 

1 

COMMON/ BLK20/ 

(5) 

321 

4 

1 

COMMON/BLK20/ 








VARIABLE NAME 

ID NUMBER 

RESTRICTION CODE 

LOCATION 

XMOLD(l) 

322 

5 

C0MMON/BLK20/ 

(2) 

323 

5 

COMMON/BLK20/ 

(3) 

324 

5 

C0MMON/BLK20/ 

(U) 

325 

5 

COMMON /BLK20/ 

(5) 

326 

5 

C0MMON/BLK20/ 

UMOLD(l) 

327 

4 

COMMON/ BLK20/ 

(2) 

32S 

4 

COMMON/BLK20/ 

(3) 

329 

4 

COMMON/ BLK20/ 

(4) 

330 

4 

COmON/BLK20/ 

(5) 

331 

4 

C0MMON/BLK20/ 

XMI(l) 

332 

1 

COMMON/ BLK20/ 

(2) 

333 

1 

COMHON/BLK20/ 

(3) 

334 

1 

C0ilMON/BLK20/ 

(4) 

335 

1 

COMMON/ BLK20/ 

(5) 

336 

1 

COMI^ON/BLK20/ 

UMI(i) 

337 

4 

COMMON/BLK20/ 

(2) 

338 

4 

COMMON /BLK20/ 

(3) 

339 

4 

COMMON/BLK20/ 

(4) 

340 

4 

COMMON/BLK20/ 

(3) 

341 

4 

COMMON/ BLK20/ 


.53 


«■ ft 4B 


Example ; 


Suppose one wishes to store on disk and later plot every 5th value of 
TM, CIA, GIB, CIC, FANG, AND RVEL 

versus TIME and that PGMODE » 3. From the preceeding table of variable 
ID numbers one obtains: 

1. TM - ID // = 21 

2. CIA - ID // = 1 

3. CIB - ID # = 2 

4. CIC - ID // = 3 

5. FANG - ID = 75 

6. RVEL - ID # * 24 

7. TIME - ID # = 79 

Since none of these variables have a restriction code number corresponding 


to FG^ff)DE * 3 the input data becomes: 

IDISK = 5 (stores every 5th value) 

IPLT = 1 (plot specified values) 

NPLTS = 6 (number of x-y plots) 

IDX • 7 (location of x variable in INDEX ) *** 

***Always place the x variable ID /‘ after the last y variable 
in INDEX i. e. IDX = NVSTOR. 

NVSTOR = 7 (Total number of variables stored NVSTOR = NPLTS + 1) 


21 

1 

Li. 


•75 

24 

79 


INDEX - 



i l■l J^ j ,i. h w 
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APPENDIX A PROGRAM LISTING 
ROUTINE PAGE 


MAIN 

NINPUT 

VIOTUT 

MINPUT 

NSETUP 

VSETUP 

MSETUP 

El'IAMOD 

I4LPMOD 

’i^LPMOD 

NETMOD 

FUNGI 

FUNC2 

INITL 

NETWHK 

NT JUMP 

SQBC 

SQMC 

SWITCH 

BCHCVP 

EPHASE 

tORPOW 

SWLOSS 

STMMF 

SANGLE 

SBGAP 

STSAN 

SERROR 

FUTURE 

MXADD 

MXMUL 

SPLPRM 

TRDIAG 

SPLINE 

SVINTP 

OUTPUT 

DVIO 

DMIO 

PDATA 

EMAPLT 

BOLPLT 

PGPLT 

YLOAD 

LINE 


A-1 

A-8 

A-13 

A-16 

A-20 

A-23 

A-27 

A-31 

A-33 

A-35 

A-37 

A-42 

A-43 

A-44 

A-48 

A-52 

A-56 

A-59 

A-62 

A-66 

A-71 

A-74 

A-75 

A-78 

A-80 

A-81 

A-82 

A-85 

A-86 

A-88 

A-90 

A-92 

A-93 

A-94 

A-95 

A-96 

A-lOl 

A-102 

A-104 

A-107 

A-110 

A-116 

A-119 

A-120 



reproducibility of the 
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INTEGRATE THE MECHANICAL POSITION LOOP FORWARD 
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